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ABSTRACT 


A  systematic  and  fundamental  study  has  been  made  of  the  microwave 
properties  of  ferroelectric  materials  with  the  perovskite  structure.  The 
quantities  singled  out  for  investigation  are  nonlinear  dielectric  constant,  micro 
wave  losses  and  electrostriction.  The  comparison  between  theory  and  experi¬ 
ment  is  discussed. 
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1.  INTRODUCTION 

A  systematic  and  fundamental  study  has  been  made  of  the  microwave 
properties  of  ferroelectric  materials  with  the  perovskite  structure.  These 
materials  have  been  chosen  due  to  their  microwave  device  potential. 1  The 
quantities  singled  out  for  investigation  were  nonlinear  dielectric  constant, 
microwave  losses  and  electrostriction.  Single  crystalline  studies  were  per- 

2 

formed  mainly  on  strontium  titanate  due  to  the  relatively  good  samples  available. 

3 

Due  to  recent  ideas  concerning  the  dynamical  suspects  of  ferroelectricity 
it  has  been  possible  to  gain  some  insight  into  the  origin  of  the  observed  micro- 

4 

wave  loss  of  materials  in  the  unpolarized  phase.  One  component  of  microwave 
loss  appears  to  result  from  the  presence  of  imperfections  in  the  crystal. 

Better  sample  preparation  should  be  able  to  minimize  this  contribution  to  the 
loss.  There  also  appears  to  be  a  contribution  to  the  loss  which  arises  from 
anharmonic  interactions  in  the  lattice.  Since  anharmonic  interactions  also 
contribute  to  the  nonlinearities,  low  loss  contribution  from  the  anharmonic 
terms  might  seem  to  imply  low  nonlinearity.  However,  since  the  anharmonic 
loss  depends  on  the  details  of  the  frequency  spectrum  and  the  nonlinearity 
depends  on  only  the  properties  of  long  waves,  low  anharmonic  loss  and  large 
anharmonic  nonlinearity  are,  strictly  speaking,  not  mutually  exclusive. 


1  "Study  of  Nonlinear  Microwave  Ferroelectric  Devices,  "  Contract  No. 

AF  33(6l6)-8022,  WADD,  Interim  Engineering  Report  No.  1  (1  March  1961  - 
31  May  1961). 

2 

Samples  were  obtained  from  National  Lead. 

3 

W.  Cochran,  in  Advances  in  Physics,  edited  by  N.  F.  Mott  (Taylor  and  Francis 
Ltd.,  London,  I960),  Vol.  9,  p.  387;  P.  W.  Anderson,  paper  given  at  All-Union 
Conference  on  Dielectrics,  Moscow,  1958  (Bull.  Acad.  Sci.  USSR,  I960),  Vol.  24; 
R.  Landauer,  H.  Juretschke,  and  P.  Sorokin  (unpublished  manuscript). . 

4B.  D.  Silverman,  Phys.  Rev.  125,  1921  (1962). 


In  section  2  of  this  report  a  procedure  is  developed  to  systematically 

interpret  the  origin  of  the  various  contributions  to  the  dielectric  constant. 

Not  only  are  anharmonic  terms  in  the  lattice  potential  energy  considered,  but 
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also  electron  deformation  terms.  One  usually  fits  the  temperature  dependence 
of  the  dielectric  constant  with  an  expression  of  the  following  form: 

C 

€  -  y-  iy~  +  Const, 
c 

In  section  2  it  is  shown  that  the  constant  cannot  be  interpreted  as  arising  solely 
from  the  presence  of  electronic  polarizibility  and  other  hard  modes  of  the  lattice; 
electron  deformation  terms  also  contribute  to  this  constant.  It  is  interesting 
to  speculate  that  large  electronic  deformations  (large  second  order  moment) 
should  contribute  to  a  large  third  order  potential  (as  Szigeti  has  shown).  Since 
the  third  order  potential  is  responsible  for  the  broadening  of  the  "ferroelectric" 
absorption  line,  the  broad  lines  that  are  observed  could  result  from  large 
electronic  deformation.  It  seems  reasonable  that  there  should  be  relatively 
large  electronic  deformation  associated  with  the  large  amplitude  of  the  soft 
mode.  Of  course  large  amplitude  also  leads  to  large  lattice  anharmonicities. 

In  section  3  the  calculation  of  the  dielectric  constant  has  been  extended  to 
include  nonlinear  effects.  Large  nonlinearity  is  essential  for  the  eventual  use 
of  ferroelectric  materials  in  many  microwave  applications. 

The  nonlinear  response  is  different  for  crystals  which  are  clamped  or 
unclamped  with  respect  to  the  probing  field.  A  measure  of  this  difference  is 
obtained  from  the  electrostrictive  constants  of  the  material.  While  there  have 
been  many  measurements  of  the  electrostrictive  and  piezoelectric  coefficients, 
there  has  been  no  theory  relating  the  results  of  measurements  to  the  microscopic 
parameters  of  the  lattice.  In  section  4  such  a  theory  has  been  developed.  It 
is  believed  that  results  of  this  calculation  will  assist  in  enhancing  our  under¬ 
standing  of  the  requirements  of  useful  materials. 


5B.  Szigeti,  Proc.  Roy.  Soc.  (London)  A252,  217  (1959). 
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2.  TEMPERATURE  DEPENDENCE  OF  THE  DIELECTRIC  CONSTANT  OF 
FERROELECTRIC  MATERIALS 
|,;i  *2.1  tetrodactlon 

Above  the  ferroelectric  transition,  the  dielectric  constant  of 
BaTiOj  and  related  isomorphs  is  large  and  exhibits  a  Curie -Weiss  type 
temperature  dependence.  It  has  been  suggested  that  this  behavior  results 
from  a  near  cancellation  of  the  short  range  restoring  forces  and  the  long 
range  driving  forces  on  the  ions. 1  Anharmonic  terms  in  the  lattice 
potential  energy  are  held  responsible  for  the  temperature  dependence  of 
the  dielectric  constant.  Since  the  Curie  temperature  is  positive*  the 
anharmonic  interactions  stabilize  the  system  above  its  transition  tem¬ 
perature.  Cochran2  has  suggested  that  the  Curie -Weiss  behavior  of  the 
dielectric  constant  results  from  the  temperature  depencehCe  of  a  long 
wave  length  transverse  optical  mode  of  the  lattice.  The  suggested 
temperature  dependence  of  this  mode  has  been  observed  by  Barker  and 

3 

Tinkham.  At  the  transition  temperature  the  frequency  of  this  mode 
should  go  to  zero  and  the  lattice  displacements  associated  with  this 
mode  become  unstable.  A  transition  is  made  to  the  ferroelectric  state. 

If  the  anharmonic  interactions  do  stabilize  the  system  above  the  transi¬ 
tion  temperature,  the  harmonic  part  of  the  frequency  associated  with 
this  mode  is  an  imaginary  quantity.  One  cannot  systematically  consider 
the  effects  of  anharmonic  interactions  on  the  equilibrium  properties  of 
a  ferroelectric  in  the  unpolarized  phase  by  expanding  the  free  energy 
in  powers  of  the  anharmonic  coupling  coefficients  about  the  harmonic 
state.  Such  an  expansion  would  lead  to  divergences.  To  avoid  this 
difficulty,  all  quantities  should  be  expanded  about  their  thermal  averages.4 
In  this  manner  one  can  systematically  derive  the  thermal  equilibrium 
properties  of  a  ferroelectric  in  its  unpolarized  phase. 

In  this  paper  we  will  derive  an  expression  for  the  temperature 
dependent  dielectric  constant  in  the  unpolarized  phase  of  a  ferroelectric 
which  manifests  a  vibrational  instability  at  a  positive  temperature.  We 
will  consider  a  Hamiltonian  such  as  Sziget^  treated  in  discussing  the 
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temperature  dependence  of  the  dielectric  constant  of  the  alkali  halides. 
Therefore,  not  only  will  third  and  fourth  order  anharmonic  terms  in  the 
potential  energy  expansion  be  considered,  but  also  election  cloud  defor¬ 
mation  through  the  presence  of  a  second  and  third  order  electronic 
moment.  It  will  be  shown  that  the  exact  linear  response  can  be  obtained 
for  this  Hamiltonian.  The  fourth  order  anharmonic  interaction  with  a 
positive  coupling  coefficient  stabilizes  the  system.  Third  order  anhar¬ 
monic  interactions  of  either  sign  raise  the  instability  temperature.  The 
largest  effect  due  to  the  presence  of  the  electronic  deformation  terms  is 
the  contribution  of  a  constant  term  to  the  dielectric  constant  at  temperatures 
high  compared  with  the  Curie  temperature.  This  indicates  that  one  cannot 
unambiguously  assign  the  origin  of  such  a  term  to  the  other  "hard"  optical 
modes  of  the  lattice. 


The  Hamiltonian  of  the  system  is 


ri  V 

-  qEVn  q°  +  y  {jgr  F(k)  +  EA(k) }  *£  +  at*  H*?w  +  a?  ) 

k  ^  L°  u>“ 

K  \J  k  k 


+  Y  fe-  G°(k)  +  E  B°(k)}  (a?  +  *°kH*°.k+*Z  > 

rg  N  uo-kkkk 

+  y  {^_  Ga(k)  +  E^Ba(k)}-~  (a?  +  aa  )<aak  +  \  > 
~  N  VfT  wa-kkk-k 


Hj  +  h2  +  \  (Pf  -  -  oE^N  q°  +  H3  +  H4  +  H5 
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The  only  difference  between  this  Hamiltonian  and  the  one  treated  by  Szigeti^ 
appears  in  the  third  term.  All  long  wavelength  transverse  optical  modes 
which  are  unstable  in  the  harmonic  approximation  have  been  collectively 
labelled  with  zero  wave  vector  and  assigned  an  imaginary  frequency.  This 
procedure  is  possible  since  the  total  number  of  such  wavelengths  is  small 
compared  with  the  total  number  of  wavelengths  of  the  optical  branch.  ®  q° 
and  p°  are  the  normal  mode  coordinate  and  conjugate  momentum  respectively 
of  this  long  wavelength  transverse  mode.  The  sums  over  k  do  not  include 

q*  o 

this  mode.  The  ak  and  ak  are  just  the  usual  creation -destruction  operators 
for  an  optical  mode  of  wave  vector  k  which  satisfy  the  following  commutation 
relations 

fak.  akil  =  [ak,  akJ  =  0  ;  [ak>  ak,]  =  *kk, 

Superscript  a  refers  to  an  acoustic  mode.  N  is  the  number  of  ion  pairs. 

The  calculation  is  performed  for  two  atoms  per  unit  cell  and  the  optical 
mode  of  this  model  is  chosen  as  the  soft  temperature  dependent  ferroelectric 
mode.  This  enables  us  to  calculate  the  temperature  dependence  of  the 
dielectric  constant  resulting  from  the  soft  mode.  Contributions  to  the 
dielectric  constant  arising  from  other  optical  vibrations  appear  just  as 
additive  constants  in  the  final  result.  The  a,  A(k)  and  B(k)  describe  the 
linear,  second,  and  third  order  moments  respectively.  The  F(k)  and  G(k) 
describe  the  third  and  fourth  order  anharmonic  douplings  respectively.  E 
is  the  macroscopic  field.  The  uk  and  wk  are  the  frequencies  of  the 
acoustical  and  optical  mode  having  wave  vector  k.  The  negative  of 
is  the  square  of  the  imaginary  harmonic  frequency  of  the  unstable  transverse 
optical  long  wavelength  mode.  This  Hamiltonian  has  been  derived  by  the 
assumption  of  a  cubic  structure  with  a  center  of  symmetry  for  the  material 
under  consideration.  The  ideal  perovskite  structure  fulfills  these  require¬ 
ments. 
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The  Hamiltonian  [Etj.  (1)J  is  transformed  by  a  unitary  trans¬ 
formation 


e'^He18  =  H  +  itH.S]  -  \  [[H.S],  S]  +  ...«**  (2) 

with  transformation  operator  S  given  by 

S  =  lJ5  I  “k  g(W  (*k+aV  (a-k  •  “k  > 

vw  k 

(3) 


-  g(k)(.a-«!k)u?k  +  <  > 

k 


This  form  has  been  chosen  for  S  so  that  H3  can  be  eliminated  from 
the  transformed  Hamiltonian  by  suitably  choosing  the  g(k). 


iHHj  +  HjJ.S]  +  H3  =  0 


(4) 


Therefore 


g(k) 


F(k)  +  EA(k) }  -j— 

*  * 


1 


(5). 


since 


[  (Hj  +  H2),S] 


.  (6) 
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Other  commutation  relations  of  use  are: 


q  tT 

[H,,SJ  ■  2i)  |^P(k)  +  EA0c)|J 

v  I™ 


TT-X  <■*&<*£ 

“k(uk  ‘ “k  1 


(7) 


°  tr  *  * 

■U1  ljp<»0  +  EA(k)|2 ^ 


) 


V'  q°  q?  q°  ♦  * 

[H4,S]  =  ”4i)  {_e_G°(k)  +  E  $.B°(k)  }  -°hg(k)(a°.  +  a?  )(aj+aa.) 

X  N  V5T  >/IT  -k  “k  ~k  -k 


(8) 


[h5,s] 


4iY  {^-  Ga(k)  +  E  —  Ba(k)  }  —  tr  g(k)  (af+a^Ha0.  +  a,°*) 
%  N  M  \Z!T  K  -k  -k  *. 


(9) 


S]  ,sj  =  2^  |i  F<k)  +EA(k)  l2--  -J  i  <a£  +  a°’)  <a°k  +  aj  f 


*  * 

0.  0  w_o  .  _o 


(10) 


-*£  I  =  *«  +  EAW  I2 - T - T  UIt,‘V(»-kt•^*,  • 

Y  jtf  K  k  k  K 


u‘(t£  -w“  ) 
k  k  k 


The  following  relations  have  been  used  to  write  the  preceding 
results  in  the  form  shown: 

S  S  £  & 

Ga(-k)  =  Ga(k)  ,  Ba(-k)  =  Ba(k) 
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With  the  use  of  these  commutation  relations,  the  transformed  Hamiltonian 
H<j.can  be  written 


n-r'lHkkk  ak  +  i)  +ItrWk(ak  ak+I)  +  i<po  "Wo  qo  }  '  oE^q! 


+  £  {  V  G°(k)  +  E  B?(k) }  -3  <a£  +  a°k) (a°k  +  a°  ) 

k  ^  Uk 

o2  o 

+  y  {  Ga(k)  +  E^-Ba(k)}-4  (»k  +  a-k)  (a-k  +  \  ) 
VN  JN  wakkk"k 


-Y  |  —  F(k)  r  EA(k)  I4 
Lj  US’ 


+  )  |  -§  F(k)  +  EA(k) 


W  ■<  > 


*  * 

/  o  .  _o  .  /  o  .  o  . 

<»k  +  i-k)(“-k+>k  > 


“k<“k  -  “k  > 


o  o  o 

+  4^  (—  G°(k)  t  E—  B°(k)}  {^2  F(k)  +  EA(k)) - ^ - 2 - r  ( a°  +  a°  111  H*) 

k  *  *  ^  I^w, 


k  k  k  k 


0  0  q°  tr  *  * 

-4Y  { Ga(k)  +  E  —  Ba(k) } {  —  F(k)  +  E A(k)} - - - - - Z“(a-k  +  ak  )(ak+a-k> 

k  *  x  .  x  (“K,1/2(< -”f> 


-i  [[H,S].§  + 


where  H"  is  H  -  (Hj  +  Hz). 
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Terms  linear  in  the  third  order  enharmonic  coupling  coefficient  F(k)  and 
second  order  moment  coefficient  A(k)  have  been  eliminated  from  the 
Hamiltonian.  Terms  bilinear  in  both  coefficients  do,  however,  appear. 

The  preceeding  transformation  has  been  performed  since  only  the  terms 
diagonal  in  the  creation -destruction  operator  that  we  have  exhibited 
explicitly  in  the  Hamiltonian  [Eq.  (11)]  will  contribute  to  the  linear 
response  of  the  system.  This  will  be  shown  in  what  follows.  Let  us 
group  the  diagonal  terms  and  label  them  Hjj.  Since  the  state  of  thermal 
equilibrium  will  be  considered,  one  can  also  set 

*k*k  *  ‘V-k  •  <»» 


Therefore 


+W(aZ**k°  + 


+ 1  (po  +  [  R  Z  ^  °°ik)  (ak  *k  +  &  *  sZ“rG“(k)  <\*£  +  i> 

k  k  k  k 


,  tt  * 

|2  /  o  o  ,  1  v 

1  - 7~  **  *k  +  F1 

k  ) 


s2lp(k) 


(13) 


2 


F(k) 


tT  a*  a  i 

—  y-o2  {<  <  + 

<  K  -4> 


n2 
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The  partition  function  of  the  system  is 


“k  '“k  k 
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H'  contains  all  terms  nondiagonal  in  the  creation -destruction  operator 
and  diagonal  terms  originating  from  £[  H'  ,  S] ,  sj  and  higher  order 
commutators. 

Expand  the  partition  sum  in  powers  of  H1 

■  r*5,,  rt,  rxA  rxs>‘lx.i.(xi'»iD„,  <k2'VHD 
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It  will  be  shown  that  in  the  limit  of  infinite  normalization  volume, 
i.  e.  N  «  ao,  Z'  will  not  contribute  to  the  linear  response.  The  free 
energy  F  can  be  written 


and  hence  the  polarization  is  given  by 


P 


_  dF  _  1  8Z 

9E  (TZ  Iff 


(16) 


(17) 


To  obtain  the  linear  response  it  is  therefore  only  necessary  to  examine 
at  most,  terms  quadratic  in  the  electric  field. 

Let  us  consider 


(18) 


The  trace  is  evaluated  in  a  representation  for  which  the  a^  are 
diagonal.  One  can  then  write 


■I  <nc 
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lnO> 
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a  a 


with  given  by  equation  J 13),  however  with  a^  a^  and 

a£  replaced  by  n£  and  i£  respectively,  n£  and  n£  are  just  the  phonal 

occupatijon  numbers  for  the  k—  acoustic  and  optical  mode  respectively. 
Note  that*  one  can  consider  the  coefficient  in  front  of  q^  as  the  square 
of  an  effective  frequency  0  which  is  dependent  upon  the  occupation 
of  all  nonzero  wave  vector  phonon  states. 


tr 
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•WK> 


It  will  be  shown  that  n£  and  can  be  rigorously  replaced  by  average 
values  n£  and  n£  to  be  defined.  If  02  (n^)  is  positive  at  all  temperatures, 
the  material  will  not  make  a  transition  to  a  ferroelectric  state  at  any 
temperature.  If  this  quantity  changes  sign  at  a  temperature  for  which 
the  phonon  occupation  numbers  can  be  replaced  by  their  classical  values, 
then 


SI2  ~  (T-Tc)  (21) 

and  above  the  transition  the  dielectric  constant  should  follow  a  Curie -Weiss 
law.  This  is  observed  for  such  materials  as  BaTiO^,  PbTiOj,  and  KNbO^. 

If  the  square  of  this  frequency  changes  sign  at  a  temperature  for  which  the 
phonon  occupation  numbers  deviate  significantly  from  the  classical  values, 
then  deviations  from  a  Curie -Weiss  behavior  of  the  dielectric  constant 
will  be  observed  before  the  material  goes  through  a  phase  transforma¬ 
tion.9  This  is  apparently  observed  in  KTa03.10 

Since  we  are  interested  in  the  linear  response,  the  partition 
function  is  expanded  to  terms  in  E2.  Terms  linear  in  E  are  dropped 
since  they  will  vanish  when  the  average  over  the  q°  coordinate  is 
performed. 
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In  the  limit  of  infinite  normalization  volume  i.e.  N.  oo,  the  result  of 
performing  the  sums  over  the  n^  is  to  just  replace  the  by  their 
thermal  average  n^  in  all  terms  except  the  first  two  appearing  in  the 
exponential,  is  given  by 


■I  v 

°k 


/ 

■piTWk^  +  j)/  2e'PlrWk(nk  +  I) 

"k 


(23) 
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This  statement  can  be  proved  as  follows*  Expand  all  terms  in  the 
partition  sum  about  the  thermal  average  n^,  except  the  first  two 
appearing  in  the  exponential.  The  first  term  in  the  expansion  is  just 


equation  (22) ,  however  with  bars  appearing  over  all  the  except 
the  first  n£  and  n£  in  the  exponential,  i.e. 


zo  =  Z  <ve 


v\ 


(24) 


R  is  easily  defined  by  comparison  of  equation  (24)  with  (22).  It  can  be 
seen  that  all  terms  quadratic  in  the  applied  field  are  proportional  to  the 
number  of  ion  pairs  in  the  normalization  volume.  In  view  of  equation  (17) 
the  polarization  is  therefore  proportional  to  the  normalization  volume 
as  expected.1^ 

The  next  term  in  the  expansion  is 

2  N  >  K> .  ,25, 

V”k  k’  iv 


All  such  terms  vanish  identically  with  use  of  equation  (23).  The  only 
nonvanishing  term  in  next  order  is 

y  <n  1.  a'R(v  -'Zj®  w 

v\  °  *'  ** 

Contributions  from  this  term  that  are  quadratic  in  the  applied  field  vary 
as  N°  and  are  therefore  neglected.  Contributions  to  the  terms  independent 
of  the  applied  field  vanish  as  N  1  when  N  •  ce.  In  a  similar  fashion,  it 
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can  be  seen  that  all  higher  order  terms  in  such  an  expansion  also 
vanish  in  the  limit  of  infinite  normalization  volume.  Therefore, 
fluctuations  about  the  thermal  values  n^  do  not  contribute  to  the 
linear  response. 

Similarly  one  can  show  that  contributions  to  the  linear  response 
from  H*  [Eq.  ( 14)]  also  vanish  in  the  limit  of  infinite  normalization 
volume.  For  that  part  of  H'  which  is  nondiagonal  in  the  creation- 
destruction  operators,  there  are  no  linear  terms  of  the  form 


<Ve)Kb 

e  1 


H*  e 


-xihd 


n> 


The  first  nonvanishing  contribution  will  come  from  a  term  of  second 
order  in  H*  ,  namely 


,  <VP>hd 

<n|e  1  U  H' 
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Contributions  from  this  term  that  are  quadratic  in  the  applied  field  vary 
as  N°  and  are  therefore  neglected.  Contributions  to  terms  independent 
of  the  applied  vanish  as  N'1  as  N*«.  In  a  similar  fashion,  it  can  be 
seen  that  all  higher  order  terms  in  such  an  expansion  do  not  contribute 
to  the  linear  response.  An  examination  of  the  terms  of  H*  that  are 
diagonal  in  the  creation -destruction  operators  will  also  show  that  they 
cannot  contribute  to  the  linear  response  of  the  system,  so  we  can  write 
[using  Eq.  (17)  and  (24)] 


Equation  (29)  is  an  exact  expression  for  the  linear  response  of  the  system 
described  by  the  Hamiltonian  given  by  equation  (1)  in  the  limit  of  infinite 
normalization  volume. 


To  evaluate  Z°  ,  matrix  elements  between  the  zero  wave 
vector  states  must  be  taken  and  then  the  sums  performed  over  all  phonon 
occupation  numbers.  Since  the  final  result  is  independent  of  the 
number  of  k  =  o  phonons,  this  procedure  can  be  simplified  by  treating 
p°  and  as  classical  variables.  One  can  therefore  write 
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Using  equation  f  29),  (31),  the  definition  of  the  dielectric  constant  • 


(<  -  1)  =  |  [v  is  the  volume  of  the  unit  cell] 


(32) 


and  the  average  phonon  occupation  numbers  [Eq.  23] 
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one  can  write 
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with 


If  the  deformation  terms  are  absent,  we  obtain 
4tt(  a^v) 


In  the  classical  limit  of  high  temperature 


02~  (T  -  Tc) 


and  the  dielectric  constant  follows  a  Curie -Weiss  law. 


(36) 


(37) 


(38) 


If  the  system  aoes  not  spontaneously  polarize  at  any  positive  temperature, 
it  is  then  stabilized  at  T=0  by  zero  point  quantum  mechanical  fluctuations. 
The  dielectric  constant  [Eq.  36]  is  independent  of  temperature  in  the 
vicinity  of  absolute  zero.  Note,  however,  that  at  intermediate  temperatures 


the  temperature  dependence  of  the  dielectric  constant  is  in  general 

9 

different  from  the  temperature  dependence  of  Barrett  ’'s  expression. 
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This  difference  arises,  since  Barrett's  calculation  is  based  upon  an 
Einstein  model  of  the  lattice  and  neglects  the  variation  of  lattice 
frequencies  with  wave  vector. 

In  the  classical  limit,  the  dielectric  constant  with  inclusion  of 
the  deformation  terms  [Eq.  34]  can  be  written  as 

«  -  1  =  AT  +  yfjfr  ■  )  [  1  +  BT  +  DT2]  .  (39) 

C 

A,B,C,D  are  temperature  independent  constants.  The  term  linear  in 
temperature  appears  since  the  presence  of  a  second  order  electronic 
moment  makes  it  possible  for  the  field  to  drive  modes  other  than  the 
soft  mode.  This  term  should  contribute  a  very  small  amount  to  the 
total  dielectric  constant.  The  presence  of  the  higher  order  electronic 
moments  also  produces  a  temperature  dependent  Curie  constant.  The 
third  order  moment  and  the  product  of  second  order  moment  with  the 
third  order  potential  contribute  to  the  linear  temperature  dependence  of 
the  Curie  constant.  The  quadratic  temperature  dependent  term  involves 
the  squares  of  the  coupling  coefficients  appearing  in  the  linear  term. 

One  can,  therefore,  conclude  that  the  largest  effect  on  the  dielectric 
constant  arising  from  the  treatment  of  the  deformation  terms  is  to 
make  the  Curie  constant  vary  linearly  with  temperature.  At  tempera¬ 
tures  high  compared  with  the  Curie  temperature,  the  dielectric  constant 
[Eq.  39]  can  be  written 

«  “  1  =  «  w —  +  F  .  (40) 

1_1C 

Therefore,  if  one  fits  the  temperature  dependence  of  the  dielectric 
constant  with  an  expression  of  the  form  of  equation(40)  over  a  wide 
range  of  temperature  above  the  Curie  temperature,  one  cannot  attribute 
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the  value  obtained  for  F  solely  to  the  presence  of  temperature  inde¬ 
pendent  infrared  active  modes.  Electronic  deformations  may  con¬ 
tribute  significantly  to  this  term. 

In  the  preceding  discussion  we  have  treated  the  crystal  as  if 

if  were  clamped,  i.  e, ,  we  have  neglected  any  effect  due  to  thermal 

expansion.  If  one  knows  the  pressure  dependence  of  the  dielectric 

constant,  compressibility,  and  volume  coefficient  of  thermal  expansion, 

the  temperature  dependence  of  the  dielectric  constant  of  the  material 

at  fixed  volume  can  be  separated  from  the  temperature  dependence 

12 

due  to  thermal  expansion.  To  obtain  an  accurate  separation  of  the 
intrinsic  temperature  effect  from  this  volume  effect,  these  three 
quantities  must  be  known  over  the  range  of  temperatures  for  which  the 
dielectric  constant  is  measured.  Such  data  is  not  available  at  present 
and,  therefore,  precludes  the  possibility  of  making  am  accurate  com¬ 
parison  of  our  results  with  experiment. 
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3.  NONLINEAR  DIELECTRIC  CONSTANT 


*  3.  l  Introduction 

The  connection  between  lattice  dynamice  and  the  dielectric  properties  of  the 

hard  ferroelectric!  (e.g.  BaTiO,  and  isomorphs)  has  been  made  both  theoretically* 
2  '  3 

and  experimentally  .  It  is  now  generally  agreed  that  the  interesting  temperature 
dependent  properties  of  ferroelectric  materials  arise  from  a  long  wavelength  trans¬ 
verse  optical  mode  which  becomes  unstable  at  the  ferroelectric  phase  transition. 
Above  the  transition  temperature  the  frequency  of  this  long  wavelength  mode  is  not 
determined  solely  by  harmonic  interactions  between  the  constituent  ions.  In  fact, 
it  appears  that  harmonic  interactions  are  not  sufficient  to  stabilize  this  mode  and 

provide  a  real  frequency  above  the  transition  temperature.  The  stabilization  of 

3 

this  mode  can  be  brought  about  by  a  consideration  of  anharmonic  interactions  . 

Since  anharmonic  interactions  perform  such  a  fundamental  role  in  providing  for 

4 

the  stability  of  the  system,  they  cannot  be  treated  in  the  usual  fashion  when  one 

considers  the  statistical  mechanical  properties  of  a  ferroelectric.  That  is  to  say, 

one  cannot  simply  expand  the  free  energy  in  powers  of  the  anharmonic  coupling 

coefficients.  It  has  been  shown  however  that  one  can  do  statistical  theory  for  a 

5 

ferroelectric  by  a  simple  modification  of  the  usual  procedure  .  In  reference  5 

it  was  shown  how  one  could  obtain,  in  the  language  of  Born-  Von  Karman  lattice 

dynamics,  the  linear  response  of  a  ferroelectric  above  its  transition  temperature. 

6 

In  the  present  paper  ,  we  consider  in  detail  the  nonlinear  response. 

In  Section  3*2  the  Hamiltonian  we  use.is  discussed.  The  Hamiltonian  is  a 

7 

slightly  modified  version  of  the  one  treated  by  Szigeti  in  discussing  the  tempera¬ 
ture  dependence  of  the  dielectric  constant  of  the  alkali  halides.  One  modification' 
consists  in  including  certain  sixth-order  anharmonic  terms  in  the  potential  energy 
expansion.  An  additional  difference  is  that  all  long  wavelength  transverse  optical 
modes  which  are  unstable  in  the  harmonic  approximation  are  collectively  labeled 
with  zero  wave  vector  and  assigned  an  imaginary  frequency.  An  expression  for 


the  temperature  dependence  of  the  nonlinear  dielectric  response  above  the  transi¬ 
tion  temperature  is  then  derived  in  Section  5.3.  It  will  be  shown  in  Section  3.  4  that 
the  temperature  dependence  of  the  nonlinear  response  we  obtain  is  almost  the 
same  as  obtained  from  the  Devonshire  phenomenological  theory.  The  Devonshire 
theory  was  quite  successful  in  describing  the  dielectric  properties  of  BaTiOj  in 
its  four  phases  with  the  use  of  one  free  energy  function.  As  suggested  by  Anderson, 
the  success  of  this  procedure  has  its  origin  in  the  nature  of  the  displacive  transition. 
If  only  a  small  number  of  modes  become  unstable  at  the  transition  it  can  easily  be 
shown  why  the  same  free  energy  function  is  valid  both  above  and  below  the  transi¬ 
tion.  This  is  discussed  in  Section  3„5.  The  deviations  from  the  Devonshire  theory 

g 

could  however  be  experimentally  significant  .  An  order  of  magnitude  estimate 

of  the  nonlinearity  and  a  comparison  with  the  measured  value  for  SrTiO^  shows 

that  it  is  reasonable  to  assume  that  enharmonic  interactions  are  responsible  for 

9 

the  observed  nonlinearity  . 

Thus  it  is  seen  that  it  is  possible  to  systematically  investigate  the  statistical 
mechanical  properties  of  a  ferroelectric  material  in  terms  of  the  Born- Von.  Karman 
lattice  dynamics  even  though  a  small  number  of  modes  have  imaginary  frequencies 
in  the  harmonic  approximation  and  that  a  lattice  dynamical  foundation  of  the 
Devonshire  phenomenological  theory  can  be  exhibited. 

3. 2  Hamiltonian 

We  shall  consider  the  following  Hamiltonian: 
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The  term*  in  the  first  line  of  Eq.  (2)  represent  the  harmonic  part  of  the  potential 
energy  due  to  lattice  displacements.  q°  is  the  normal  coordinate  of  that  optical 
mode  designated  by  zero  wavenumber  and  which  has  a  purely  imaginary  frequency 
associated  with  it.  The  specification  of  k  =  0  means  that  we  consider  a  mode 
with  wavelength  long  compared  with  the  lattice  parameter  but  small  compared 
with  the  sample  size.  There  are  a  large  number  of  such  modes*®.  However, "they 

are  only  a  very  small  fraction  of  the  total  number  of  modes  and  are  therefore 

o  a 

treated  as  a  single  mode.q  and  q  are  the  normal  coordinates  of  wave  number  k 

K  K 

for  the  optical  and  acoustic  mode,  respectively.  The  second  and  third  lines  of 
Eq.  (2)  represent  the  third-and  fourth-order  anharmonic  contributions  to  the 
potential  energy,  respectively.  N  is  the  number  of  ion  pairs.  The  calculation 
is  performed  for  two  atoms  per  unit  cell  and  the  optical  mode  of  this  model  is 
chosen  as*the  soft  temperature  dependent  ferroelectric  mode.  The  second  term 
on  the  right  of  Eq.  (1)  represents  the  first-order  effect  of  the  electronic  dis¬ 
tortions  plus  the  effect  of  the  displacement  of  the  ions  as  a  whole.  E  is  the 
macroscopic  electric  field.  The  first  and  second  lines  of  Eq.  (3)  are  the  second* 
and  third-order  effect  of  electronic  deformations,  respectively.  Equation  (4) 
represents  the  sixth-order  anharmonic  contribution  to  the  potential  energy.  The 
remaining  undefined  quantities  are  constants.  The  summations  over  k  are  under¬ 
stood  to  be  for  all  values  of  k  different  from  zero.  We  have  not  explicitly  in¬ 
cluded  in  the  Hamiltonian  terms  arising  from  the  kinetic  energy  since  this  will 
not  contribute  to  the  dielectric  response.  In  this  same  vain,  we  shall  neglect 
in  the  calculations  that  follow  all  multiplicative  constants  that  do  not  effect  the 
final  result.  The  Hamiltonian  has  been  simplified  by  the  assumption  of  a  cubic 
structure  and  a  center  of  symmetry  for  the  material  under  consideration.  These 
assumptions  concerning  the  structure  are  consistent  with  the  perovskite  struc¬ 
ture  in  the  unpolarized  phase.  All  other  terms  in  the  Hamiltonian  which  have 
been  neglected  can  be  shown  to  make  contributions  smaller  than  the  order  of 
approximation  to  which  we  eventually  work. 


26 


3.3  Dielectric  Response  Abort  Transition  Twagntart 


The  work  of  the  preaant  section  is  divided  into  two  parts.  In  Part  A,  we 

neglect  the  sixth-order  enharmonic  contributions  to  the  Hamiltonian.  This 

5 

Hamiltonian  then  corresponds  exactly  with  the  one  used  previously  in  obtaining 

an  expression  for  the  linear  response.  We  exactly  calculate  the  first  two  terms 

2 

in  a  power  series  in  E  which  represents  the  dielectric  response.  At  temperatures 

well  above  the  transition  temperature  the  sixth-order  anhar moni citie s  can  be 

shown  to  contribute  relatively  little  to  the  nonlinear  response**  (their  major  effect 

is  to  introduce  a  term  proportional  to  E*  in  the  dielectric  response).  This  then 

is  the  reason  we  neglect  them  in  this  first  part.  In  Part  B  we  neglect  the  effect 

of  electron  deformations.  The  sixth-order  anharmonicities  are  treated  in  a  less 

2 

exact  way.  The  first  three  terms  in  a  power  series  in  E  for  the  dielectric 
response  are  obtained.  Below  the  transition  the  sixth-order  anharmonicities 
are  needed  to  stabilize  the  system.  In  a  latter  section  we  discuss  the  behavior 
below  the  transition.  The  results  of  this  section  will  then  be  of  great  importance. 

A.  >  0 

Our  Hamiltonian  consists  of  the  first  three  terms  of  Eq.  (1): 
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=  X 


-  e  EN*  q 


+  <XL 


(5) 


The  partition  function  for  this  system  is  given  by: 
oo  so  oo 

za  =  r  "i?  i  ^  $  £  d%  •**  <  •  >  v  •  “> 

-00  .00  .00 

o  a 

where  (3  =  1/k^T.  The  integrals  over  dq^  and  dq^  may  be  readily  performed. 

In  doing  this,  we  make  use  of  the  following  symmetry  properties  of  the  coefficients 
in  the  Hamiltonian: 
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Equation  (6)  then  becomes 
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Above  the  transition  temperature,  in  the  absence  of  the  term  quartic  in  q^,  the 
system  described  by  Eq.  (7)  would  be  unstable  if  the  enharmonic  terms  in  the 
potential  energy  were  neglected.  That  these  terms  do  indeed  stabilise  the  system 
is  made  evident  by  rewriting  Eq.  (7)  in  the  following  form: 
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o2 

*0 


(xk  +  *1 


(11) 


/  0  f  1  _  o2  o2  0^-1  o4|\ 

#xp  \"pLtwo  %  -•  EN  q0  +  N  FqoJj  * 

where  the  "  renormalized"  frequency  u®  is  proportional  to  J T-T  and  is  thus 
o 

real  for  T  >  T  .  «  is  given  by: 

C  w 


“o  =  ”“o  +  Np  Z  (x,k+  **k  "  vk)  ’  <12) 

k 

Tc  is  a  positive  number  (the  Curie  temperature)  and  can  be  obtained  from  Eq. 
(12).  It  is  shown  later  [  Eq.  (17)]  that  the  field  independent  parts  of  the  product 
expression  of  Eq.  (11  )  which  are  proportional  to  q°2  cancel.  The  N  dependence 
of  the  remaining  terms  is  such  that  they  can  in  no  way  affect  the  stabilization  of 
the  system. 

Since  we  wish  to  obtain  an  expression  for  the  dielectric  response  (or  the 
polarisation)  as  a  power  series  in  E,  we  evaluate  the  expression  for  the  partition 
as  a  power  series  in  E.  E  is  found  in  both  the  exponent  and  the  function  T^  of 
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Eq.  (11 ).  The  expansion  of  the  former  is  Itraightforward.  The  expansion  of 
1-4 

T^  *  is  also  trivial.  One  must  however  then  expand  a  function  of  the  form 
l"I{l  +  Aj'E  +  . . .  }.  This  can  be  accomplished  by  making  use  of  the  following 
result,  the  proof  of  which  is  easily  given: 

+\E  +  BkE2  +  ckE3+DkE4+...> 

= 1  +  e(Zv  +  E2{IBk  +  Z  W> 

k  k  k  <k* 

+  E'  <Z  S  +  Z  <Vk. +  V  V + Z  W  V  > 

k  k<k'  k<k'  <k" 

+  E<t  Z  Dk+ Z  <Akck.  tAk.ck+BkV  i 

k  k<k' 

+  Y  (AkAk' Bk»  +  Ak"Bk  +  W'  ) 

k<k'  <k" 

+  Y  AkAk'  Ak"Ak"»  ^  + 

k<k»  <k"<k" ' 

■  1  +  PE  +  QE2  +  RE3  +  SE4  +  . . . 

Upon  making  the  necessary  expansions,  Eq.  (11 )  then  becomes: 
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where  3C^  »  GJq^ qQ2  +  N  *  F q^4 .  The  quantise e  P,  Q,  R,  S  are  defined 

in  terms  of  the  A^,  B^,  C^,  by  Eq.  (7).  The  A^,  B^,  C^,  are  given 
below: 
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where 
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We  have  neglected  the  terms  of  odd  power  in  E  in  Eq.  (14)  since  these  can  be 

o  o 

shown  to  be  of  odd  power  in  q^  and  hence  vanish  when  the  integral  over  q^  is 

evaluated.  We  have  also  changed  the  lower  integration  limit  in  Eq.  (14)  and 
continue  to  neglect  unimportant  multiplicative  factors  which  make  no  contribution 
to  the  dielectric  response.  Before  performing  the  integral  in  Eq.  (14)  a  word 
should  be  said  concerning  the  dependence  upon  N  of  our  results  (hence  an  explana¬ 
tion  of  the  terms  written  as  0  ( N  ^q^)  of  Eq.  (15)).  The  final  quantity  (the 
polarization)  we  will  calculate  must  be  proportional  to  N,  which  is  a  very  large 
number.  Thus  in  what  follows,  we  neglect  any  contributions  which  are  propor¬ 
tional  to  N  *  as  compared  to  unity.  Therefore  we  show  explicitly  only  those 
terms  which  contribute  to  the  desired  N-  dependence  of  the  final  result.  One 
final  step  must  be  taken  before  we  evaluate  the  integral.  We  must  simplify  the 
product  expression  in  the  integrand  of  Eq.  (14).  Bearing  in  mind  our  previous 
discussion  about  powers  of  N,  this  can  be  accomplished  by  use  of  a  result  equiv¬ 
alent  to  Eq.  (13)  [  N  1  replaces  E  as  the  small  parameter]  .  When  this  is 
carried  out  we  find: 
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Substituting  this  and  tbs  valuss  of  P,Q,R,S  givts  for  Eq.  (14): 
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where  the  and  are  constants  (independent  of  q^  and  N)  and  can  be  obtained 

by  directly  comparing  Eq.  (18)  with  Eq.  (14)  and  making  use  of  Eqa.  (13).  (15) 

0 

and  (lq).  The  integrals  over  dq^  are  now  readily  performed.  This  is  accom¬ 
plished  by  making  use  of  the  following  result: 


=  a"n  T(n  +  ~  )  [  1  -  ya‘2(n  +  ~-)(n +“ )  +  0  (y2)] 
when  yo  <  <  1  . 

Thus  we  obtain  the  following  expression  for  Z  : 
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The  dielectric  response  is  now  obtained  from  Eq.  (20)  by  means  of  the  fol¬ 
lowing  definition: 
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which  gives 
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Instead  of  the  definition  of  the  dielectric  response  given  by  Eq.  (22)  we  could 
have  introducted  the  so-called  differential  dielectric  response  by 
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This  does  not  change  the  first  term  on  the  right  of  Eq.  (23)  but  would  place  a 
factor  of  three  in  front  of  the  second  term.  In  this  section  we  continue  to  use  the 
definition  of  the  dielectric  response  given  by  Eq.  (22). 

That  part  of  the  dielectric  response  independent  of  E  will  be  referred  to 
as  the  linear  dielectric  response.  All  other  terms,  which  depend  on  E,  will  be 
called  the  non-linear  contribution  to  the  dielectric  response.  The  linear  dielectric 
response  for  the  Hamiltonian  under  consideration  here,  has  been  discussed  in 
detail  in  reference  5.  The  result  obtained  here  corresponds  to  the  classical  limit 
of  their  result  [  their  Eq.  (39)  ]  •  From  Eq.  (23)  we  directly  obtain,  neglecting 
all  terms  down  by  a  factor  of  N  ^  compared  to  unity: 
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where  by  F^(y)  we  mean  that  part  of  y  proportional  to  x  and  where 
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We  now  write  down  the  expression  for  that  part  of  the  non-linear  response 
2 

proportional  to  E  .  Note  that  since  c  -1  should  not  be  proportional  to  N.  one  must 

worry  about  the  factor  of  N  in  front  of  the  second  term  on  the  right  of  Eq.  (23). 

.  2 

However,  it  is  found  that  the  terms  independent  of  N  in  b  -  f  a  cancel  exactly 
and  that  the  lead  contribution  goes  as  N  *  ,  cancelling  the  factor  of  N  in  the 
appropriate  term  of  Eq.  (23).  Thus  we  again  obtain  the  proper  dependence  on 
N  and  our  procedure  is  valid. 
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and  where  the  1^,11^,111^  IV^  are  given  by  Eq.  (16). 

If  electronic  deformations  are  neglected,  the  only  coefficients  In  Eq.  (27) 
that  do  not  vanish  are  E^,  Ej  and  F^.  Equation  (27)  for  the  non-linear  response 
then  becomes 
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From  Eq.  (27)  it  appears  that  the  major  effect  of  including  the  electronic  deforma¬ 
tions  is  to  add  a  term  E^F^  to  H'  .  Then  H!  in  Eq.  (29)  is  replaced  by  H  where 
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The  Hamiltonian  we  coneider  in  this  part  is: 

kb  =  Ki  -  EN^  +Ky  (32) 

i 

We  treat  3C^  approximately  0C^  is  sufficient  to  stabilize  the  system)  and  calculate 
only  the  lowest  order  effect  of  this  term.  In  doing  this  we  shall  repeatedly  make 
use  of  the  following  results: 


<F:  rk'  °k>  ~  <ak2> 

(G;  Pk;  Qk;  S fc;  Tk)  -  (F2;  Fk2;  G^) 


(33) 


By  this  we  mean  that  any  of  the  quantities  on  the  left  are  of  the  same  order  of  mag¬ 
nitude  as  any  of  the  quantities  on  the  right.  That  is  to  say,  the  fourth-order  potential 
coefficients  are  the  same  order  of  magnitude  as  the  square  of  the  third*x>rder  poten¬ 
tial  coefficients  etc.  In  this  same  vain  we  could  also  write: 
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In  what  follows  then,  we  neglect  any  terms  whose  magnitude  is  smaller  than  the 
sixth-order  potential.  Thus  we  will  neglect  terms  which  are  proportional  to  GF  etc. 


The  partition  function  for  the  system  is  given  by  Eq.  (6)  with  X. 
by  3Cg.  Denoting  by  X^'  that  part  of  3Cj  which  directly  involves  the  q 
3C  3'  3  3Cj  "  N"2Gq^)  we  simplify  the  calculation  by  writing 
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The  integrations  over  dq  and  dq  in  the  partition  function,  Z  ,  may  then  be 

k  k  B 

performed  in  a  straightforward  manner.  After  this  is  done,  again  neglecting  all 
unnecessary  multiplicative  factors,  we  findt 


C  „  °  f af  1  02  02 

J  dt>o  **>“  qo  - 


4  o  -1  o2  -2 

e*  EN*  qQ  +  N  rqQ  +  N  Gq 


t{sM 


,  .  a  /  -1  -2  o2  .  ,  -2-1  o4 

1  +  P  (N  p  qp  +  N  p  qp  v] 


Jl)  - 


?]) 


(35) 


where 


(34) 


(37) 

pk  s  Vwk2 :  V’  V^k2 

We  now  proceed  as  inlfert  A.  To  explicitly  show  the  stabilisation  of  the  system 

above  the  transition  temperature  by  the  higher  order  potential  energy  terms,  we 

o2 

remove  from  the  product  terms  those  contributions  which  depend  on  qQ  and  thus 
"renormalize"  thefrequancy  ,  To  within  the  approximation  which  we  work, 

Eq.  (35)  is  rewritten  as  below: 


“k  ■  -<38k  +  J%+V 


vk = 


-  K +  pk> 


and 


_  /  ad 

V“k  1 


>:VS 


o2 


42 


/ 


'B 


•  j*  **°9  exP  f  P  [j^2 

-00  ' 


=o2  o2  4  o  -1  o4  , 

e*  EN  qQ  +  N  FqQ  + 


n 


k  ,,tp  ^n_1 


qQ2  (Xk  +  \  ’  Vk} 


•]}?{* 


(38) 


where  the  "renormalized"  frequency  Mq  is  given  by 
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Note  that  the  renormalized  frequency  is  proportional  to  -Tc  +  T  +  irT^  (v  «  1). 
Both  product  terms  of  Eq.  (38)  are  simplified  as  described  previously.  When 
this  is  done  the  second  line  of  Eq.  (38)  is  replaced  by  {l  +  o  N‘*  q°4}  where 


a  = 


(40) 


and  r  is  defined  by  Eq.  (17).  The  partition  function  is  finally  written  in  the  follow¬ 
ing  form  consistant  with  the  previous  approximation: 
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The  reason  for  writing  Z  in  this  form  will  be  made  evident  in  Section  31.5. 


The  partition  function  as  given  by  Eq.  (41)  may  now  be  directly  evaluated. 

We  again  evaluate  the  partition  function  as  a  power  series  in  even  powers  of  E  (the 
odd  powers  vanish  identically).  Since  the  renormalized  frequency  is  positive,  we 
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do  not  need  the  higher  power  terms  in  q^  (i.  e. ,  ~  q^,  q^)  in  the  exponent  to 
stabilise  the  system  so  that  in  the  evaluation  of  the  integral  we  may  treat  them 
as  perturbations.  Z  is  then  found  to  be  of  the  same  form  as  Eq.  (20).  However, 
we  shall  now  calculate  one  additional  term  in  this  series.  Therefore  a  term  NJcE 
is  added  to  the  right  of  Eq.  (20).  The  dielectric  response  is  now  obtained  from 
Eqs.  (22)  and  (23),  the  latter  being  modified  by  the  addition  of  a  term 
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Thus  the  dielectric  response  is  written  as 


c  -  1  =  J  +  KE2  +  LE4  +  . . . 
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with 


J  =  C/(T  -  Tc  +  S-T  ) 


2.4 


K  =  (M  +  H"  T)/(T  -  Tc  +  vT  ) 


(43) 


L  =  U/(T  -  Tc  +  eT2)6  +  V/(T  -  Tc  +  <rT2)7 


and  where 


=  -2kB  [h  2  \]  /[if  2  (\  +  '‘k  -  Vk>] 


H"  =  4  a 


4  4 

k  e*  T 
B  c 


08 


(44) 


44 


u 


-6G 


V 


48F 


o  14 
"0 


(44) 


C,  Tc  and  M  are  defined  in  Eqs.  (26)  and  (30).  This  result  is  exact  to  terms  of 
order  N~*,  that  is,  terms  which  vary  as  N"*  compared  to  unity  are  neglected. 
The  point  is  again  made  that  this  result  neglects  contributions  which  involve  any 
coefficient  whose  magnitude  is  smaller  than  the  sixth-order  enharmonic  coupling 
coefficients. 


3.4  Comparison  with  the  Devonshire  Theory 

In  the  previous  section  we  developed  expressions  for  P  (or  equivalently 

12 

«  -  1)  as  a  series  in  powers  of  E.  We  now  perform  a  simple  inversion  of  these 
series,  that  is,  we  express  E  as  a  series  in  powers  of  P  so  that  the  preceding 
results  may  be  directly  compared  with  the  results  of  Devonshire^  phenonenological 
theory.  The  discussion  in  this  section  is  restricted  to  temperatures  above  the 
transition  temperature. 


Consider  first  the  results  of  Section  3.3-A.  Combining  Eq.  (25)  and  Eq.  (29) 
with  the  modification  as  indicated  by  Eq.  (31),  the  polarisation  may  be  written: 
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where  the  first  term  on  the  right  of  Eq.  (45)  has  been  simplified  by  neglecting  the 

small  term  linear  in  the  temperature  and  assuming  T»  T  so  that  T/T  -  T  *  1. 

c  c 
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can  now  be  assumed  to  consist  of  three  parts;  firstly,  a  contribution  from  the 
electronic  polarizability  and  from  the  other  hard  modes  which  until  now  had  been 
neglected  and  lastly  *  a  contribution  BC  arising  from  the  effect  of  electronic 
deformations.  In  Eq.  (45)  we  have  written  p  for  N  *F.  Inverting  Eq.  (45)  by 
a  self  consistent  technique  results  in  the  following: 
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Equation  (46)  can  be  cast  into  the  following  form  (again  assuming  T»Tc): 
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where 
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Let  us  now  consider  the  results  of  Section  3 >3«b.  It  follows  from  Eq.  (43) 

that 
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Combining  the  results  of  Eqs.  (47)  and  (49)  we  conclude,  that,  to  the  order  of 
approximation  to  which  we  work  the  result  of  considering  the  total  Hamiltonian 
as  given  by  Eq.  (1)  is  an  expression  of  the  following  form: 
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where 


(T  =  <r  +  7 

:  _  (5D 

H  =  H  +  H"  -  H'  . 

A  factor  of  H’  has  been  subtracted  from  the  right  hand  side  of  the  equation  defining 
H  so  that  the  contribution  to  H  arising  from  the  term  proportional  to  r  is  not 
counted  twice. 

Thus  both  <r  and  H  involve  contributions  from  the  electronic  polarizability,  o:.' 
the  hard  modes,  the  electronic  deformations  and  the  sixth-order  potential  energy 
terms.  He  also  has  a  contribution  which  comes  from  the  square  of  the  fourth-order 
potential  (i.  e. ,  ~  r  as  given  by  Eq.  (17).  We  remark  at  this  point  that  this  result 
implies  that  in  considering  the  linear  response,  if  one  writes  the  dielectric  con¬ 
stant  as 


c 


(52) 


that  « j  involves  contributions  from  the  electronic  polarizability,  the  hard  modes, 
the  electronic  deformations  and  from  the  sixth-order  potential  energy  terms.  This 
point  was  considered  previously  in  reference  5. 


The  expression  for  E  vs  P  above  the  transition  temperature  can  now  be 
directly  compared  with  the  work  of  Devonshire.  Devonshire  writes  an  expression 
of  the  form  of  Eq.  (50)  but  with  .  r  =  H  =  0.  Eq.  (50)  indicates  that  the  result 
of  the  most  general  calculation  would  be  of  the  form: 


That  the  terms  involving  ir  and  H  are  not  negligible  and  may  indeed  be  experimentally 

observed  will  be  discussed  in  Section  3.  6.  Corrections  to  the  Devonshire  theory  as 

8 

proposed  here  have  previously  been  considered  by  Triebwasser  and  by  Kozlovs 


ski?3. 
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3.  5  Dielectric  Response  Below  the  Transition  Temperature 


Up  to  this  point,  we  have  assumed  that  the  temperature  was  above 
transition  temperature.  It  is  known,  however,  that  even  below  the  transition 
temperature  an  expression  of  the  form  given  by  Eq.  (50)  or  Eq.  (53)  can  be 
used  to  explain  the  observed  dielectric  response.  We  now  wish  to  briefly 
indicate  how  this  follows  directly  from  the  lattice  dynamics  of  the  problem. 

In  what  follows  we  shall  neglect  the  effect  of  electronic  deformations  so  that 
we  again  consider  the  Hamiltonian  as  given  by  Eq.  (32).  We  work  to  the  same 
degree  of  approximation  as  discussed  previously.  We  further  assume  that, 
although  this  Hamiltonian  was  derived  under  the  assumption  of  certain  sym¬ 
metry  conditions  (Section  3.2)  which  are  no  longer  strictly  valid  below  the 
transition,  the  form  of  the  Hamiltonian  is  the  same.* 


-q2 

When  T<  T  ,  the  square  of  the  renormalized  frequency  is  no 

longer  positive  and  the  system  again  appears  to  be  unstable.  Further,  terms 
- 1  o4 

of  the  form  N  Fq„  (F>0)  are  also  unable  to  stabilize  the  system  since  due 

0  3 

to  electrostrictive  effects  F-*F'<  0.  Thus  it  is  left  to  terms  of  the  form 
N-2GqQ^(G>0)  arising  from  the  sixth-order  anharmonic  potential  energy  to 
stabilize  the  system.  Thus  the  thermodynamic  treatment  of  this  section 
must  differ  fundamentally  from  that  of  the  previous  sections. 


The  reason  for  rewriting  the  partition  function  Zg  in  the  form  given 
by  Eq.  (41)  will  now  be  made  evident.  We  may  look  upon  the  factor  multiplying 
-p  in  the  exponential  of  Eq.  (41)  as  an  effective  Hamiltonian  (the  system  being 
stabilized  by  the  term  N_2Gq°^).  Below  the  transition  temperature  the  polar¬ 
ization  is  directly  related  to  the  displacements  q^  so  that  we  may  again  simply 
obtain  an  expression  for  E  vs  P.  From  Eq.  (41)  we  have: 


Keii  =  2  '  e*  EN^0  +  N-^V-ap-1)  +  N_2q°6G  +.  .  .  .  (54) 
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If  the  system  is  in  equilibrixm  the  free  energy  must  be  a  minimum. 
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eff 


=  0  . 


(55) 


In  essence,  we  evaluate  the  partition  function  by  means  of  the  method  of  steepest 
descents.  Equation  (55)  then  defines  the  saddle  point.  Substituting  Eq.  (54)  in 
Eq.  (55)  then  gives 
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It  can  be  simply  shown  that 
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so  that  we  have 


E  =  p 


|u°2  e*  +  p3  j^4e*  ^  (®*  "  QP  +  P*  ^6e*  •  •  • 


(58) 


which  is  identical  to  Eq.  (50).  It  thus  appears  that  the  basic  reason  that  one  free 
energy  function  works  both  above  and  below  the  transition  temperature  is  intimately 
related  to  the  fact  that  we  have  considered  only  a  single  unstable  mode  and  that  the 
other  modes  remain  unchanged. 


3.  6  Conclusion 

It  has  been  shown  how  one  can  obtain  an  expression  for  the  nonlinear  response 
of  a  ferroelectric  in  terms  of  the  coupling  parameters  of  the  Born- Von:  Karman 
lattice  dynamical  theory.  Up  to  the  present  all  other  calculations  along  these  lines 
have  been  essentially  single  particle  calculations  and  have  hence  neglected  certain 
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contribution*  made  manifest  by  the  lattice  dynamical  treatment.  The  role  of  elec¬ 
tron  cloud  deformation*  has  also  been  investigated  and  shown  to  contribute  to  the 
temperature  dependence  of  the  nonlinear  response.  It  has  also  been  shown  how  the 
assumption  of  one  unstable  mode  used  in  conjunction  with  the  same  Hamiltonian  both 
above  and  below  the  transition  is  equivalent  to  Devonshire's  assumption  that  one 
phenomenological  free  energy  function  suffices  to  describe  the  behavior  both  above 
and  below  the  transition.  Only  order  of  magnitude  estimates  are  possible  at  present 
to  justify  the  origin  of  the  observed  effects.  The  detailed  calculation  that  has  been 
presented  is  of  interest  however,  since  the  quantities  involved  also  have  implications 
for  the  observed  optical  behavior  of  ferroelectric  materials.  Indeed,  as  Szigetti. 
has  pointed  out,  widths  of  absorption  bands,  the  presence  of  sidebands,  etc,  are 
determined  by  essentially  the  same  coupling  parameters  that  are  involved  in  the 
temperature  dependence  of  the  dielectric  response. 
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4.  ELECTROS TICT1VE  CONSTANTS  OF  FERROELECTRIC  MATERIALS 


4. 1  Introduction 

The  piezoelectric  coefficients  have  been  related  to  the  microscopic  prop¬ 
erties  of  the  crystalline  lattice  by  Bom  and  Huang  in  their  classic  text  on  lattice 
dynamics.*  This  theory  is  applicable  when  the  induced  stress  in  a  material  is 
proportional  to  the  applied  field.  If  however,  every  atom  in  the  lattice  is  at  a  cen¬ 
ter  of  symmetry,  this  linear  effect  is  absent  and  the  stress  is  proportional  to  the 
square  of  the  field,  hi  this  paper  it  will  be  shown  how  one  can  relate  this  quadratic 
behavior  to  the  microscopic  properties  of  the  lattice.  The  quadratic  dependence  of 
stress  on  field  has  been  called  electrostriction  by  a  number  of  investigators  and  this 
is  the  terminology  that  will  be  adopted  here. 

Electrostrictive  effects  are  generally  small  and  difficult  to  measure  in  most 
materials.  Ferroelectrics  however  exhibit  pronounced  electromechanical  behav¬ 
ior.  These  materials  undergo  a  transition  to  a  polarized  phase  below  a  certain 
temperature.  Above  the  transition  temperature  they  are  paraelectric,  with  a 
large  dielectric  constant  which  obeys  a  curve  Weiss  law.  One  may  classify  these 
materials  as  belonging  to  one  of  two  categories;2  those  that  are  piezoelectric  in 
the  unpolarized  phase  and  those  that  are  electrostrictive  in  the  unpolarized  phase. 

The  piezoelectric  ferroelectrics  have  been  recently  discussed  by  Cochran  and 
will  not  concern  us  in  the  present  work.  Barium  titanate  and  its  isomorphs  be¬ 
long  to  the  second  category  since  in  their  high  temperature  cubic  phase  each  ion 
is  at  a  center  of  symmetry.  The  electrostrictive  properties  of  these  materials 
are  of  interest  for  several  reasons.  It  appears  that  the  explanation  of  the  ob¬ 
served  first  order  phase  transition  in  BaTiO,  requires  the  presence  of  suitable 

3  4 
electrostrictive  terms  in  the  expression  for  the  free  energy  of  the  material. 

Electrostrictive  behavior  also  leads  to  the  difference  between  the  value  of  the 

5 

clamped  and  free  nonlinear  dielectric  constant  in  the  paraelectric  phase.  One 

4 

pertinent  assumption  of  Devonshire's  theory,  is  that  the  piezoelectric  coeffi¬ 
cients  obtained  for  the  material  in  the  polarized  phase  are  related  simply  to  the 
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electrostrictive  coefficients  of  the  material  in  the  unpolarized  phase.  The  micro- 

6 

scopic  basis  for  such  assumption  has  been  discussed  by  Anderson.  In  this  sense 
one  could  therefore  say  that  the  electromechanical  behavior  of  barium  titanate  in 
all  phases  is  electrostrictive. 

The  calculation  relating  the  electrostrictive  coefficients  to  the  microscopic 

coupling  parameters  of  the  lattice  also  presents  some  new  and  interesting  problems 

not  encountered  in  the  linear  theory.  Toupin  has  stated  that  in  general  one  cannot 

consider  the  stress  tensor  as  a  polynomial  in  the  electric  field  and  elements  of 

7 

infinitesimal  strain  measure.  He  has  demonstrated  that  this  assumption  violates 
the  invariance  of  the  stored  energy  to  rigid  rotations.  To  the  order  we  will  work, 
however,  it  will  not  be  necessary  to  augment  the  strain  tensor  with  an  antisymmetric 
part  representing  rotations  in  the  medium.  We  can  forego  this  complication  since 
we  treat  a  lattice  in  which  each  atom  is  at  a  center  of  symmetry.  It  should  there¬ 
fore  be  emphasized  that  the  expressions  obtained  for  the  electrostrictive  coeffi¬ 
cients  in  this  paper  are  not  valid  for  the  quadratic  response  of  a  material  which 
also  exhibits  a  linear  response.  Toupin  also  raises  the  following  interesting  point. 

If  the  electric  field  and  polarization  are  not  parallel  to  each  other,  the  Maxwell 
stress  tensor  in  the  dielectric  medium  is  not  symmetric.  Hence  if  the  medium 
is  to  be  in  static  equilibriumthere  must  be  an  additional  stress  system  whose, 
antisymmetric  part  is  equal  and  opposite  to  the  antisymmetric  part  of  the  Max¬ 
well  stress  tensor.  One  can  call  these  additional  stresses,  local  stresses.  In 
the  calculation  presented  in  Section  4. 2  it  Will  be  shown  that  Maxwell  stresses  and 
local  stresses  can  be  identified  and  that  both  are  in  static  equilibrium  with  each 
other. 

It  might  appear  at  first  that  in  order  to  calculate  an  effect  that  is  quadratic 
in  the  macroscopic  field  one  would  need  to  work  within  the  framework  of  the  theory 
of  finite  elasticity.  This,  however  is  not  necessary  since  even  though  the  macro¬ 
scopic  field  depends  on  the  ionic  displacements,  it  also  depends  on  the  geometry 


54 


6f  the  sample  through  the  depolarisation  factor  and  therefore  the  field  and  die- 
placement  are  independent  quantities.  One  can  therefore  work  to  any  given  order 
in  each  independently. 

4.  2  Electro  strictive  Constants 

Expressions  for  the  electro  strictive  constants  are  obtained  after  a  compari¬ 
son  has  been  made  between  the  macroscopic  and  microscopic  theories.  In  this 
section  we  will  therefore  consider  the  elastic  equation  of  motion  for  an  electro- 
strictive  medium  and  the  microscopic  lattice  equations  of  motion.  These  equations 
will  be  put  into  a  form  which  will  facilitate  a  comparison  between  them.  The 
elastic  equation  of  motion  is 


(1) 


p  is  the  mass  density  of  the  medium,  u  is  the  a  component  of  the  displacement 

th  a  th 

in  the  medium.  S  is  the  av  component  of  the  stress  tensor.  X  is  the  v 
o»  T  v 

component  of  a  general  point  in  the  medium. 

The  constitutive  relation  between  stress,  strain  and  field  is 
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6A  *6A  + 


2 


d 
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E  E 

6A  6  A 


(2) 


We  assume  that  the  piezoelectric  effect  is  absent  since  each  ion  is  considered  tdc 
be  at  a  center  of  symmetry.  The  C  .  .  and  d  .  .  are  the  elastic  constants  and 
electrostrictive  constants  respectively.  Eg  is  the  6^  component  of  the  macro¬ 
scopic  field,  s.  is  the  6A^  component  of  the  strain  tensor  which  is  related  to 
the  displacement  gradients  as  follows. 
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(3) 
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Substituting  Eq.  (2)  and  Eq.  (3)  into  Eq.  (1)  we  obtain 
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The  symmetry  of  the  elastic  and  electro strictive  constants  have  been  used  to  ob¬ 
tain  this  result.  The  displacements  and  macroscopic  field  are  now  Fourier 
analyzed. 


u 

a 


2iri  y  •  r 
e 


(5) 


Er$*E<,y  •2"iy' r 


(6) 


If  Eq.  (5)  and  (6)  are  substituted  into  Eq.  (4)  and  the  usual  manipulations  per¬ 
formed  to  obtain  uncoupled  motion  in  the  harmonic  approximation,  one  obtains 

P  V  =  4’2  X  Cav  ,  64  yv  y4  "«y 
V,  6,  A 


+ **2  d.v . 64 JI w  n4<?'  -*>  <7» 

v,  6,  A 


e 
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6(y)is  the  Dirac  delta  function.  It  ie  seen  from  Eq.  (7)  that  modes  of  wave  vector 
yare  no  longer  uncoupled.  Different  modes  are  coupled  due  to  the  presence  of  the 
electrostrictive  term.  We  will  obtain  expressions  for  the  eiectrostrictive  con¬ 
stants  from  a  comparison  of  Eq.  (7)  with  the  microscopic  equations  of  motion. 
Therefore  we  turn  to  an  examination  of  the  microscopic  equations  of  motion.  In 
what  follows,  the  rigid  ion  model  is  used  for  simplicity.  The  generalisation  to 
more  sophisticated  models  within  the  adiabatic  approximation  should  be  straightforward. 

The  lattice  equation  of  motion  in  the  notation  of  Born  and  Huang  is 


se  i 

"VW  * 


I 


Vkk'1  V£,  + V. 


f'.K.p 


-n 


6  (<  i' f«)  u  (i')tt  (i") 

Topy  'kkf  k«'  p'k* '  y'k» 1 


(8) 


f.k'.p 

i-.k-.y 


m,  is  the  mass  of  the  k**1  type  atom,  u  (*)  is  the  a**1  component  of  the  displace- 
ment  of  the  k  type  atom  in  the  i  cell.  The  $  .  and  ^  are  second  and  third 


opy 


th 


order  coupling  coefficients  respectively,  q.  is  the  ionic  charge  of  the  k  type 
til  ^ 

atom,  is  the  a  component  of  the  macroscopic  field.  The  harmonic  inter¬ 
action  has  been  split  up  into  a  short  and  long  range  part  as  discussed  by  Born 
and  Huang.  The  macroscopic  field  «  therefore  includes  contributions  from  the 
long  range  interaction  between  ions  and  $  is  determined  by  the  short  range  in- 
ter  actions  between  ions. 
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After  application  of 


w  (*) 
a  k 


1 

<s/N 


^  wa(k|y)  e2ffi  y  **k* 

7 


with 


(9) 


(10) 


the  equation  of  motion  Eq.  (8)  can  be  written  in  the  following  form. 
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Vk' 


|y)  + 


€  _ 
ay 


(ID 
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where  we  have  set 


Vk£> 


(nW 
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6  I**') 

vap'k  k' 1 
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<b  V 

^apy  k  k'  k"' 


(13) 
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<U) 


y  ij> 


A(y'  +"y’'-  y)  =  1 
a  0 

One  method  of  identifying  the  electrostrictive  coefficienta  in  Eq.  (7)  with  the  micro* 
acopic  parameter  a  of  Eq.  (11)  ia  the  following.  Let  ua  aaaume  that  up  to  time  zero, 
the  nonlinear  interactiona  in  both  of  theae  equations  are  not  preaent.  At  time  zero 
they  are  inatantaneously  awitched  on.  The  displacementa  are  aet  equal  to  un¬ 
perturbed  diaplacementa  plua  small  connections  due  to  the  preaence  of  the  non¬ 
linear  interactions.  We  will  treat  Eq.  (11)  first.  Set 


for  y*  +  y"  -  y  =  0  or  a  reciprocal  lattice  vector, 
all  other  values  of  yr ,  y*',  and  y 


(15) 


w  (k|y)  *  w  (kjy)  +  v  (k|y)  (16) 

a  oc  ci 

Subscript  zero  refers  to  the  displacements  in  the  absence  of  the  third  order  an- 

harmonic  coupling  parameter.  For  an  infinitesimal  time  after  the  enharmonic 

8 

interaction  has  been  turned  on,  we  may  write. 


Vjkly)  .  -  l  e2'**  '  1  <”  &  Vk'  &  + 

ro,  V » P 


€  — 
ay 


(17) 
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y 1  y"  mk'p 
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A(y'  +y"-y)  (18) 
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u  _ 
ay 


The  elaetic  Eq.  (7)  can  be  treated  in  a  similar  fashion.  The  displacement 
is  set  equal  to  an  unperturbed  part  plus  a  small  correction 


u  _  =  u  _  +  v  _ 
ay  oay  ay 


In  the  same  spirit  in  which  Eq.  (17)  and  (18)  have  been  written  we  write 


(19) 


p  u  _  =  -4ir 
oay 
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C  e .  y  y  .  u._ 
a»,tA  »  A  6y 


(20) 


v,  6,  A 


p  v  _ = 4ni  y 
ay  Li 


av ,  6A 


d7  dy"  Egyff  Ea-,  y^  6(7<  +  7  -  y)  (21) 


v,  6,  A 


Equations  (17)  and  (20)  are  just  the  equations  that  Born  and  Huang  use  to  obtain 
expressions  for  the  elastic  constants.  Had  we  included  the  piezoelectric  terms 
in  Eq.  (20);  these  equations  could  be  used  to  determine  expressions  for  the  piezo¬ 
electric  coefficients  also.  To  obtain  expressions  for  the  electrostrictive  coefficients 
we  must  compare  Eq.  (18)  with  (21).  Before  this  comparison  is  possible  however, 
Eq.  (18)  must  be  written  in  the  limit  of  long  waves;  the  quadratic  dependence  on 
the  macroscopic  field  must  be  exhibited;  the  cepter  of  mass  motion  for  the  unit 
cell  must  be  obtained.  In  what  follows  we  describe  in  some  detail  how  this  is  done. 

To  obtain  the  motion  of  the  center  of  mass  of  the  unit  cell  we  multiply  both 
sides  of  Eq.  (18)  by  Vm^,  sum  over  the  contents  of  the  unit  cell  and  divide  by  the 
sum  of  the  masses  in  the  unit  cell. 
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I  V. 
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2iriiy‘-x(S^)+y,,-5fi.v)l 


A(y»  +  y*  -  y)  (22) 


For  wavelength*  large  compared  with  the  atomic  interaction  diatance,  the 
exponential  in  Eq.  (22)  is  expanded. 


.*<?•<£> ♦y-sfr&i  _  1  +  », ,  * ... 


The  first  term  in  this  expansion  appearing  on  the  right  of  Eq*  (22)  is 


X  2  •^v(k”^vk'|7',v‘!'|?,4(?+?-^  <24> 

Z,*”*  y,y"mk'P 
k  nkf-y 

k 

and  this  term  vanishes  due  to  the  following  translational  invariance  condition 


V  .  ,o  m  n  .  n 

Z  ♦«0y{kk'k">  "° 
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Contributions  to  Eq.  (22)  from  the  next  term  in  the  exp&naion  do  not  vanish, 
and  we  therefore  write 


v,<7)=-~i—  1  X  ^  Daf>Y(£  ”  E->Vk^>wOY<k«|y*>(?.S<™> 
Lm k  y.y”  mk'P 


nk -  y 
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+  yB-x(J„J)]  A(7*+r-y)  (26) 


If  the  following  interchange*  are  performed  in  the  term  arising  from  the  second 
quantity  in  the  square  brackets 

y  ^y11,  m^n,  k' s*kM,  P  s=sy 

we  can  write 


..  -  2rri 
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Let  us  next  decompose  the  displacements  into  the  normal  modes 


«'oa(kfr»*  2<«(klj7|Q(j'» 


(28) 


y 

The  equation  of  motion  for  the  Q(j )  are 


Q(J)  +  «2(J)  txj)  =  ^  «>!]■> 
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We  can  write 
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K.  ** 

y 

where  q(j )  ie  the  eolution  of 


(30) 


qf)+«2(J)q(jy)  =  0 


(31) 


Substituting  (28)  and  (30)  into  (27)  we  obtain 
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(32) 


Finally, the  Kronecher  delta  function  goes  over  to  a  Dirac  delta  function  in  the 
limit  of  infinite  periodicity  volume. 

6<r>  1  F~Nv«  J  ** 

a  7 
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I 


We  obtain 
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where  use  hee  been  made  of 


v  *  >mv  v  (y) 
ay  a  a  7 


(34) 


E  »  ‘/R’v  E  (y) 

ay  a  a 

▼  is  the  volume  of  the  unit  cell. 

a 

Equation  (33)  can  now  be  compared  directly  with  Eq.  (21)  and  the  following 
expression  obtained  for  the  electrostructive  coefficients 


2  2  bf-  J  - 

j',j"  midp  km  J  3 


nk"  y  k"' 
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/  mk  _  ,om. 
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These  coefficients  are  subject  to  the  following  symmetry  relations 

d  ad  ad 

av,  6A  va,  4 A  av,  A6 

and  in  general  there  are  thirty-six  independent  constants. 


(35) 


(36) 
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At  we  have  mentioned  in  the  introduction,  Toupin  has  pointed  out  the  fol¬ 
lowing  inter  eating  fact.  If  the  electric  field  and  polarisation  in  a  medium  are 
not  parallel  to  each  other,  the  Maxwell  etreee  tensor  in  the  dielectric  medium 
ie  not  symmetric.  Hence  if  the  medium  ie  to  be  in  static  equilibrium  there 
must  be  an  additional  etreee  system  whoee  antisymmetric  part  ie  equal  and  op¬ 
posite  to  the  antisymmetric  part  of  the  Maxwell  etreee  teneor.  The  previous 
calculation  should  include  the  Maxwell  stresses  since  we  have  been  working  to 
terms  second  order  in  the  field.  The  Maxwell  stress  tensor  ie  bilinear  in  the 
electric  field  and  polarization.  Since  the  induced  polarization  is  proportional  to 
the  field  however,  the  Maxwell  stress  teneor  is  essentially  quadratic  in  the  electric 
field  variables.  The  symmetry  of  the  electrostrictive  coefficients  Eq.  (36)  how¬ 
ever  always  leads  to  symmetric  stresses  in  the  medium.  This  guarantees  the 
rotational  invariance  of  the  energy  function  of  the  medium*  It  is  therefore  ap¬ 
parent  that  in  the  calculation  just  described.  Maxwell  stresses  and  local  stresses 
have  been  lumped  together.  This  will  be  shown  in  what  follows. 

The  third  order  coupling  parameter  appearing  in  Eq.  (18)  can  be  written  as 
the  sum  of  two  terms. 


_  o  mn 
Dapv(k  k'k") 


_N  .o  m  n  »  c 
DQpY<k  k'  kn)  +  D, 


o  m  n 

aPy  k  k'  k"' 


(37) 


Superscript  N  indicates  the  short  range  non  coulomb  contribution  to  the  third 
order  coupling  parameter.  This  term  can  include  the  effects  of  covalency,  and 
general  n  body  interactions.  Superscript  c  refers  to  two  body  coulomb  interactions 
and  therefore  k,  k',  and  kn  must  each  be  equal  to  one  of  two  types  of  ions  for 
Dapy(k  k*  k") to  be  non-vanishing.  Let  us  restrict  attention  to  the  coulomb  term, 
and  write  the  equation  of  motion  Eq.  (18)  for  this  term. 
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The  two  terms  appearing  in  Eq.  (38)  can  be  combined  if  the  following  relations 
are  used  to  rewrite  the  first  term. 
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Using  Eq.  (13)  we  can  write 
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Since  we  want  to  calculate  the  motion  of  the  center  of  mass  of  the  unit  cell,  one 
must  multiply  both  sides  of  Eq.  (40)  by  earn  over  k  and  then  divide  by  the 

sum  of  the  masses  in  the  unit  cell.  We  obtain  the  following  equation 
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after  substitution  of  Eq.  (42)  into  Eq,  (40)  and  the  dm  of 
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The  dot  over  the  equality  sign  in  Eq.  (44)  indicates  that  we  have  only  exhibited  the 
macroscopic  field  term  on  the  right  hand  side  of  this  equation.  The  inner  field 
term  which  we  neglect,  may  be  grouped  with  the  short  range  forces.  If  we  per¬ 
form  the  sum  over  p,  k'  and  k  in  Eq.  (43)  and  let  the  periodicity  volume  become 
infinite  we  obtain 


pVay  »  2ni  dy"  6(y»  +y*  -  y)  +... 


This  term  on  the  right  of  Eq.  (45)  is  just  the  Maxwell  stress  term  that  appears  in 
the  equation  of  motion.  One  can  show  this  in  the  following  way.  The  elastic 
equation  of  motion  can  be  written 
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Components  of  the  Maxwell  stress  tensor  are  given  by 


M 

S  «  E  P 
ay  ^Y 


Using  the  following  Fourier  decomposition  of  electric  field  and  polarisation 
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it  can  be  shown  that  the  Maxwell  stress  term  Exhibited  e;:plicitly  in  Eq.  (46)  is 
identical  with  the  term  in  Eq.  (45).  Hence  we  have  shown  that  the  microscopic 
lattice  equations  of  motion  that  have  been  used  to  obtain  expressions  for  the 
electrostrictive  constants,  contain  the  Maxwell  stress  terms. 
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4.  3  Electrostriction  of  Ferroelectric  Materials 


The  piesoelectric  and  electrostrictive  behavior  of  ferroelectric  materials 

10 

has  been  treated  extensively  in  the  past  Many  of  the  measurements  have  been 
made  on  materials  in  the  polarised  phase.  In  this  phase  the  materials  are  piezo¬ 
electric  exhibiting  a  linear  relation  between  stress  and  applied  field.  The  inter¬ 
pretation  of  measurements  in  the  polarised  phase  is  complicated  by  the  presence  of 
ferroelectric  domains.  In  the  unpolarized  or  cubic  phase  the  material  is  electro¬ 
strictive,  i.  e. ,  there  is  a  quadratic  dependence  of  the  stress  on  electric  field. 

4 

Devonshire's  treatment  of  the  electromechanical  properties  of  BaTiO^  assumes 
that  the  piezoelectric  coefficients  in  the  polarized  phase  are  related  simply  to  the 
electrostrictive  coefficients  in  the  unpolarized  phase.  Verification  of  Devonshire's 
theory  has  formed  the  basis  for  much  of  the  work  on  electrostriction.  To  the 
extent  that  the  Devonshire  theory  is  valid,,  the  theory  presented  in  this  paper 
should  form  a  basis  for  understanding  the  microscopic  origin  of  the  piezoelectric 
•ad  electrostrictive  behavior  o i  B.TiOj  and  related  i.omorph..  In  any  .vent,  the 
microscopic  theory  presented  in  Section  4. 2  should  be  valid  for  application  in  the 
unpolarised  phase  of  these  materials. 

11  12 

Measurements  on  single  crystal  barium  titanate  and  strontium  titanate 
have  been  made  above  the  Curie  temperature  as  a  function  of  temperature.  Over 
the  range  of  measurement,  these  materials  are  paraelectric  having  a  large  diel¬ 
ectric  constant  which  obeys  a  Curie-Weiss  law.  The  measured  temperature 
variation  of  the  electrostrictive  coefficients  d,  is  given  by 


d  **  (T  -  Tc)2  (493 

T  is  the  Curie  temperature  obtained  from  high  temperature  dielectric  measure- 

c 

mentf .  This  temperature  variation  is  in  agreement  with  the  result  predicted  from 
Devonshire's  phenomenological  theory.  That  this  temperature  variation  also 
results  from  the  microscopic  expressions  obtained  is  shown  in  the  following  way. 
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The  frequency  of  the  low  lying  "ferroelectric"  mode  of  SrTiO^  has  an  observed 
temperature  dependence  given  by1* 


. o  .  .  1/2 

«(,.)*  (T  -Tc) 


(50) 


Such  a  temperature  dependent  mode  should  also  be  present  in  BaTiO,  and  its 

6  14  ^ 

isomorphs  *  .  This  is  the  mode  which  dominates  the  low  frequency  dielectric 

behavior.  This  mode  should  also  dominate  the  electrostrictive  behavior  due  to  its 

relatively  low  frequency  and  large  oscillator  strength1*.  Therefore,  in  the 

expression  for  the  electrostrictive  coefficients,  Eq.  (35),  contributions  are 

neglected  from  all  optical  modes  other  than  the  soft  ferroelectric  mode.  The 

result  can  be  written 
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The  temperature  dependence  of  the  soft  mode  frequency  (Eq.  50)  therefore  leads 
to  the  observed  temperature  dependence  of  the  electrostrictive  coefficients. 

Since  we  have  the  electrostrictive  coefficients  in  terms  of  the  microscopic 

parameters  of  die  lattice,  let  us  see  whether  these  parameters  can  be  reasonably 

chosen  to  reproduce  the  measured  value.  The  room  temperature  value  of  d  , 

13  11  11 

for  SrTiO^  as  reported  by  Winter  and  Rupprecht  will  be  calculated*  This  value  is 
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This  coefficient  ie  written 
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This  expression  can  be  evaluated  simply  with  use  of  the  linear  chain  model  of  a 
16 

ferroelectric.  For  the  chain  model,  we  have 


« (1)  =  -  ;  q.  =  -  ne/b 

*/2 


«(2)  =  —  ;  q,  »  ■  ne/b 

VT  2 


n  is  the  number  of  electronic  charges  on  the  ions,  e  is  the  electronic  charge, 
b  is  a  constant  which  includes  effects  of  electronic  polarizability  and  makes  the 
and  q  essentially,  effective  chargee.  Two  unit  cells  of  the  chain  model  are  shown 
in  Figure  1.  After  the  sums  have  been  performed  in  Eq.  (53)  we  obtain  the 
relatively  simple  expression  for  d^,  jj. 


"ini 

*n  w  (  f)  v^ 

a  is  a  third  order  coupling  coefficient,  d  is  a  lattice  parameter  shown  in  Figure  1 
Values  have  been  substituted  for  all  of  these  quantities  (except  n)  in  a  manner 
discussed  previously.  The  value  of  n  was  then  adjusted  so  that  the  value  of  (55) 
was  equal  to  the  observed  vtflue  (52).  In  this  fashion  one  obtains  the  reasonable 
value  of  n  =  2  electronic  charges,  which  is  further  substantiation  for  the  treatment 
presented. 
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4*4  4pp*°dix, 

In  tibia  section  we  sketch  how  the  symmetry  relations  (Eq.  36)  are  obtained*. 
The  symmetry  of  the  electrostrictive  coefficient  (Eq.  35)  upon  the  following 
interchange 

6  **  A  (56) 

is  obtained  simply  by  noting  the  following  interchanges 


k"  S=s  k"  ;  k'ss  k«  ;  p  s*  y  ;  )  ^  j'  ;  m  *  n  (57) 

and  by  using  the  translational  invariance  condition  (Eq.  ’25). 

The  symmetry  of  die  electrostrictive  coefficient  upon  the  interchange 

o  5=4  v  (58) 

is  somewhat  more  involved.  Upon  interchange  of  a  and  v  in  equation  (35)  we  obtain 
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The  rotational  invariance  condition  connecting  third  order  and  second  order 
coupling  parameters  can  be  written  ^ . 
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Using  equations  (35),  (59),  and  (60),  and  the  defining  equations  (12)  and  (13),  one 
can  write 


dva,  &  "  dav,  6A 


^  ^  {DoY&)apv  +  Dpo(k'k->i 


yv  "  Dvv(£E")apa 


nk"v  kf 


DPv4S«,5Y^  x 


</(°)  ^ 


w9 


(61) 


If  the  contribution  to  the  macroscopic  field  is  separated  from  the  contributions 
linear  field  and  noncoulomb  forces  in  the  second  order  coupling  parameter. 
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It  can  be  shown  that  the  term  on  the  right  of  equation  (61)  vanishes  thereby 
proving  the  validity  of  equation  (36). 
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5.  OPTICAL  VIBRATIONS  IN  PBROVSKITE  MATERIALS 


It  is  now  generally  agreed  that  the  interesting  paraelectric  properties  of 
the  perovskites  result  from  the  temperature  dependence  of  a  low  lying  transverse 
optical  mode.  Indeed,  one  can  understand  qualitatively  the  origin  of  the  Curie - 
Weiss  behavior  of  the  dielectric  constant,  nonlinear  response,  electromechanical 
behavior,  and  microwave  loss  mechanisms  in  terms  of  a  general  lattice  dynamical 
approach.  A  more  detailed  and  quantitative  understanding  of  these  and  other 
phenomena  is  frustrated,  however,  by  our  lack  of  knowledge  concerning  the  actual 
ionic  motions  associated  with  the  soft  ferroelectric  mode  and  other  optical  modes 
of  the  material.  There  have  been  various  attempts  to  interpret  optical  data  and 
assign  motions  to  the  optically  active  vibrations;  however,  none  of  the  work  to 
date  is  conclusive.  In  this  we  briefly  describe  an  approach  to  this  problem  which 
has  yielded  results  of  some  interest. 

The  perovskite  structure  has  five  atoms  in  the  unit  cell  and  therefore  a 
total  of  fifteen  branches  in  its  spectrum.  Three  of  these  are  of  the  acoustic 
type.  Three  of  these  are  of  a  torsional  type  and  optically  inactive.  Of  the  nine 
remaining,  six  are  transverse  (doubly  degenerate)  optically  active  vibrations 
and  three  are  optically  inactive  longitudinal  vibrations.  The  long  wavelength 
optically  active  modes  can  be  characterized  by  the  relative  motions  of  the  five 
inequivalent  sublattices.  Due  to  the  cubic  symmetry  of  the  structure  and  the 
fact  that  three  of  the  sublattices  are  composed  of  oxygen  atoms  it  is  easily  seen 
that  if  one  restricts  attention  to  vibrations  along  one  of  the  cartesian  axes, 
there  are  in  general  seven  sublattice  couplings.  Since  there  are  three  observed 
optical  frequencies  (SrTiOj),  1  one  can  choose  three  of  the  sub  lattice  couplings 
different  from  zero  and  have  a  completely  specified  problem.  One  can  then 
solve  for  the  sublattice  couplings  and  also  for  the  relative  motion  of  the  sub¬ 
lattices  for  each  mode.  This  calculation  has  been  carried  out  for  the  sublattices 
as  labelled  in  Fig.  2,  and  for  the  following  couplings; 


*W.  G.  Spitzer,  R.  C.  Miller,  D.  A.  Kleinman  and  L.  E.  Howarth,  Phys.  Rev. 
126,  1710  (1962). 
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FIGURE  2 


kjs  02  -  Ti  coupling 

k2*  02  -  03  and  02  -  04  coupling 

k3:  Sr  -  03  and  Sr  -  04  coupling 

The  equations  of  motion  are 


ml*l 

* 

-  2kjXj  +  2kjX2 

m2*2 

* 

-  2kjX2  +  2kjXj  -  8k2x2  +  4k2x 

m3S3 

S 

-  4k3x3  +  *k3Xg  -  41£2X3  +  4^2X 

m4*4 

= 

-  4k3x4  +  4k3Xg  -  4k2x4  +  4k2x 

“5*5 

s 

-  8k3x6  +  4kjXj  +  4k3x4 

The  solutions  of  these  equations  are  presently  being  analyzed;  however,  one 

interesting  feature  that  automatically  comes  out  of  the  calculation  is  the  following. 

The  high  frequency  mode  is  a  vibration  of  the  02  sublattice  against  the  03  and  04 

sublattices  which  are  rigidly  fixed.  The  two  lower  frequency  motions  essentially 

involve  a  fairly  rigid  motion  of  the  entire  oxygen  framework  (02,  O y  and  04 

sublattices).  This  is  consistent  with  the  observed  ionic  displacements  from  the 

cubic  phase  of  ferroelectric  materials  of  the  perovskite  structure  as  the  material 
2 

becomes  polarized.  The  oxygen  framework  remains  relatively  rigid.  Since 
this  indicates  that  oxygen- oxygen  couplings  are  relatively  strong  it  is  therefore 
reasonable  that  the  highest  frequency  mode  is  a  vibration  between  the  oxygen 
sublattices.  Since  the  two  other  modes  involve  little  relative  motion  of  the  oxygen 


2G.  Shirane,  R.  Pepinsky,  and  B.  C.  Frazer,  Phys.  Rev.  97,  1179  (1955); 
B.  C.  Frazer,  H»  Donner,  R.  Pepinsky,  Phys.  Rev.  100,  ~T45  (1955); 

H.  T.  Evans,  Mass.  Inst.  Technol.  Lab.  Insulation  Research  Tech.  Rept. 
No.  58  (1953). 
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sublattices,  it  is  again  reasonable  that  their  frequencies  are  both  relatively 

distant  from  the  highest  frequency  as  has  been  observed.  That  this  high  frequency 

is  a  common  feature  of  the  perovskite  structure  can  be  seen  by  the  reflectivity 

3 

spectra  illustrated  in  Fig.  3.  through  16. 

Similar  calculations  have  been  performed  for  a  four  sublattice  model  in 
which  bond  angle  couplings  have  been  introduced  between  the  titanium  and  03 
and  O.  oxygen  atoms,  The  observed  torsional  infrared  inactive  frequency  was 

’  4 

used  to  completely  determine  the  four  couplings.  Again  similar  results  were 
obtained  for  the  high  frequency  mode. 


3 

The  reflectivity  data  presented  in  Fig.  3  has  been  obtained  by  G.  Rupprecht 
prior  to  his  leaving  Raytheon  Company, 

4 

We  are  indebted  to  R.  A.  Cowley  for  giving  us  the  value  of  this  infrared  inactive 
frequency  which  he  obtained  from  inelastic  neutron  scattering  measurements. 
See  Phys.  Rev.  Letters  9,  159  (1962). 
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FIGURE  3 
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FIGURE  5 
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FIGURE  6 
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FIGURE  9 
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FIGURE  10 


FIGURE  II 
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FIGURE  12 
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FIGURE  13 
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FIGURE  14 


FIGURE  15 


FIGURE  16 


6.  ELASTIC  CONSTANTS  OF  KTaQ, 

The  elastic  stiffness  constant,  Cn,  of  KTa03  was  measured  as  a  function 
of  temperature  in  order  to  determine  whether  KTaOj  exhibits  a  phase  trans¬ 
formation  similar  to  that  observed  in  SrTi03  at  110*K.  Since  the  microwave 
losses  increase  rapidly  at  such  a  transformation  its  presence  might  seriously 
limit  the  device  applications  of  a  material. 

The  elastic  constant  was  determined  by  measuring  the  velocity  of  propaga¬ 
tion  of  a  longitudinal  ultrasonic  pulse  along  a  [100]  crystallographic  axis  of  a 
single  crystal  since  vjong  a  V^l i/p*  apparatus  and  the  techniques  used 
are  identical  to  those  described  by  Bell  and  Rupprecht  for  SrTiQ3*  so  they  will 
not  be  discussed  here.  Measurements  have  been  made  from  room  temperature 
to  liquid  nitrogen  and  no  transition  has  been  observed  in  this  range.  The  results 
are  shown  in  Fig.  17.  A  change  in  ultrasonic  attenuation  which  was  observed  was 
not  reproduced  in  other  samples. 


IR.  O.  Bell  and  G.  Rupprecht,  Phys,  Rev.  129,  90  (1963). 
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Temperature  dependence  of  C|| 
in  KTa03  from  -200°  C  to  25#  C 


Temperature  (%) 
FIGURE  17 


7 .  LOW  TEMPERATURE  DIELECTRIC  CONSTANT  OF  SrTiO? 

The  dielectric  constant  of  single  crystal  SrTi03  has  been  measured  from 
2*K  to  60  *K,  These  measurements  were  made  on  a  small  disc  at  a  frequency 
of  100  Kc/sec.  The  results  are  shown  in  Fig.  18  and  are  in  substantial  agree¬ 
ment  with  the  values  reported  by  Weaver.  1 2  The  data  were  fitted  with  an  ex¬ 
pression  of  the  form  derived  by  Barrett  for  the  low  temperature  dielectric 
constant  of  perovskite  type  crystals.  The  best  fit  was  obtained  with  the 
expression 


e 


_ C _ 

T./2  +  coth  T./2T  -  T 
l  1  c 


with 

C  =*  8. 72X104 

Tc  ■  43*K 

T  j  =  1 10*K 


1H.  E.  Weaver,  J.  Phys.  Chem,  Solids  H_,  274  (1959). 

2J.  H.  Barrett,  Phys.  Rev.  86,  118  (1952). 
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8.  LOW  TEMPERATURE  DIELECTRIC  MEASUREMENTS  OF  KTaQ^ 

The  dielectric  constant  and  loss  tangent  of  single  crystal  KT&O3  were 
measured  from  2*K  to  80*K.  The  capacitance  and  equivalent  parallel  resistance 
of  a  parallel  plate  capacitor  of  circular  cross  section  with  a  KTa03  dielectric 
were  measured  at  100  Kc/sec  in  order  to  determine  the  dielectric  constant  and 
loss  tangent.  The  results  are  shown  in  Fig.  19  and  Fig.  20  .  The  dielectric 
constant  data  are  fitted  with  Barrett's  expression 


Tj/2  coth  Tl/2T  -  Tc 


with 


C  *  6.025  X  104  T  *  13. 4*K  T,  3  42*K 

c  1 

The  agreement  with  K-band  data  at  the  higher  temperatures  is  quite  good. 

The  large  peak  in  the  loss  tangent  in  the  vicinity  of  Tc  indicates  that  there 
may  be  an  actual  ferroelectric  transition  at  this  temperature.  A  comment  should 
be  made  on  the  magnitude  of  the  loss  tangent  at  the  higher  temperatures.  At  25*K 
the  equivalent  series  resistance  of  the  capacitor,  its  contacts  and  leads  is  of 
the  order  of  one  ohm  or  less.  A  large  portion  of  this  resistance  is  probably 
due  to  the  lead  and  contact  resistance  (the  lead  resistance  was  measured  to  be 
about  0. 1  ohms  but  the  resistance  of  the  silver  electrodes  on  the  sample  has  not 
been  measured).  Thus,  the  actual  losses  of  the  dielectric  at  temperatures 
above  about  10*K  are  much  smaller  than  measured  here.  These  errors  are  of 
1.0  consequence  in  the  vicinity  of  the  peak  in  the  loss  tangent  since  the  series 
resistance  is  as  large  as  several  hundred  thousand  ohms  and  the  lead  and  contact 
resistance  remain  essentially  constant  with  temperature.  Hence  the  peak  should 
actually  be  even  more  pronounced  than  the  one  shown.  For  this  reason  no 
attempt  was  made  to  fit  the  experimental  curve. 
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Temperature  (°k) 
FIGURE  19 


Loss  tangent  of  KTaO 


Temperature  (*k) 
FIGURE  20 
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